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Abstract 

We observed a linearly sliced area of the Cygnus Loop from the north-east to the 
south-west with Suzaku in seven pointings. After dividing the entire fields of view 
(FOV) into 119 cells, we extracted spectra from all of the cells and performed spectral 
analysis for them. We then applied both one- and two-component non-equilibrium 
ionization (NEI) models for all of the spectra, finding that almost all were significantly 
better fitted by the two-component NEI model rather than the one-component NEI 
model. Judging from the abundances, the high-A;T e component must be the ejecta 
component, while the \ow-kT e component comes from the swept-up matter. Therefore, 
the ejecta turn out to be distributed inside a large area (at least our FOV) of the 
Cygnus Loop. We divided the entire FOV into northern and southern parts, and 
found that the ejecta distributions were asymmetric to the geometric center: the 
ejecta of Si, S, and Fe seem to be distributed more in the south than in the north 
of the Cygnus Loop by a factor of ~2. The degree of ejecta-asymmetry is consistent 
with that expected by recent supernova explosion models. 

Key words: ISM: abundances - ISM: individual (Cygnus Loop) - ISM: supernova 
remnants - X-rays: ISM 
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1. Introduction 

The asymmetry of a supernova (SN) explosion is considered to be a key point to under- 
stand the mechanism driving an SN explosion. The spectropolarimetry of young SNe revealed 
a non-spherical core of the SN explosion. The degree of departure from spherical symmetry 
is considered to be anti-correlated with the mass of the hydrogen envelope, i.e., the remain- 
ing hydrogen shell of the progenitor star (e.g., Leonard et al. 2006). In the extreme case of 
departure from spherical symmetry, for example Type-IIn SN 1998S, the major-to-minor axis 
ratio of the explosion core is expected to be larger than 2.5, suggesting a highly asymmetric 
explosion mechanism, possibly mediated by jets (Wang et al. 2001). From theoretical aspects, 
recent SN explosion models generally show aspherical ejecta distributions may also be produced 
by convection (e.g., Herant et al. 1994; Burrows et al. 1995; Janka & Muller 1996) as well as 
hydrodynamical instabilities (e.g., Kifonidis et al. 2006; Burrows et al. 2007; Janka et al. 2007) 
during SN explosions. 

After a SN explosion, the ejecta initially expand without deceleration (so-called "free 
expansion") behind the forward shock wave, which sweeps up the interstellar medium (ISM). 
When the mass of the swept-up matter becomes comparable to that of the ejecta, the expansion 
is decelerated by the swept-up matter. The deceleration initiates a reverse shock propagating 
into the ejecta. Both the swept-up matter and the ejecta are hot enough to emit X-rays. 
The swept-up matter forms a shell structure, while the ejecta fill its interior. Therefore, we 
can attempt to address the ejecta distributions by observing the structure of supernova rem- 
nants (SNRs) in X-rays. In fact, ejecta-dominated X-ray emission are detected in many SNRs 
(e.g., Cas A: Hughes et al. 2000; Tycho: Decourchelle et al. 2001; Vela: Tsunemi et al. 1999; 
G292.0+1.8: Park et al. 2001). Recent deep Chandra X-ray observations of Cas A SNR have 
revealed a bipolar structure of Si-rich ejecta (Hwang et al. 2004; Laming et al. 2003), suggesting 
an asymmetric explosion with bipolar jets. 

The Cygnus Loop is a nearby (540 pc: Blair et al. 2005) proto-typical middle-aged 
(~10000 yrs) SNR. A number of pieces of evidence that the SN explosion of the Cygnus Loop 
occurred in a preexisting cavity (e.g., McCray 1979; Levenson et al. 1999), suggesting that the 
Cygnus Loop was produced by a core-collapse SN. Since it is an evolved SNR, the ejecta-material 
is embedded in swept-up matter. Miyata et al. (1994) observed the northeast (NE) shell of the 
Loop with ASCA, and revealed the metal deficiency there. Recent Suzaku observations of 
this region confirm the metal deficiency there (Miyata et al. 2007). Although the reason for 
the metal deficiency is not yet understood, it can be safely stated that contamination of the 
ejecta-material is negligible in the rim regions. In contrast to the rim regions, ASCA detected 
evidence of Si-, S-, and Fe-rich plasma at the center portion of the Cygnus Loop (Miyata et 
al. 1998), which is thought to be ejecta. The relative abundances of the ejecta support that 
the idea that the Cygnus Loop was the result from a core-collapse SN, and that the progenitor 
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Fig. 1. ROSAT HRI image of the entire Cygnus Loop. The Suzaku FOV (P8, P12, P13, P14, P15, P16, 
and P17) are shown as white rectangles. 

mass is estimated to be 25 M Q based on the central portion of the Cygnus Loop (Miyata et 
al. 1998). Recent XMM- Newton observations across the Cygnus Loop from the NE rim to 
the SW rim revealed ejecta distributed inside of ~ 0.85 R s of the Cygnus Loop, where R s is 
the shock radius (Tsunemi et al. 2007). The relative abundances inferred for the total ejecta 
are almost consistent with those expected for the core-collapse SN, whose progenitor mass 
is 13 M (Katsuda & Tsunemi 2008), which can be considered to be more reliable than that 
estimated from the ASCA observation, since the coverage of the total FOV of the XMM- Newton 
observations is about an order of magnitude larger than that of the ASCA observation, giving 
a more representative sample of the ejecta abundances to compare with theoretical models. 

Since the Cygnus Loop is very large in apparent size (2.°5x3.°5: Levenson et al. 1997; 
Aschenbach & Leahy 1999), it is an ideal target to study the spatial distribution of the ejecta 
material. In order to reveal ejecta structures in this SNR, we observed the Cygnus Loop from 
the NE rim to the SW rim with Suzaku (Mitsuda et al. 2007). We selected our fields of view 
(FOV) so as to cover just southern regions of the XMM-Newton observation path since the 
distribution of Si-, S-, and Fe-rich ejecta is suggested to be elongated toward the south of the 
Loop (Tsunemi et al. 2007). 
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Table 1. Information of observations of the Cygnus Loop and Lockman Hole. 



UDS. 1JJ 


Coordinate (RA, DEC) 


Roll 


Ubs. Date 


Effective Exposure 




Cygnus Loop 






501028010 (P8) 


314.005, 31.464 


58.°5 


2006.05.13 


4.7 ks 


ouiuzyuiu (±lZ) 


313.751, 31.263 


58.7 


zUUo.5.l)9 


16. Z KS 


501030010 (P 13) 


313.498, 31.061 


58.9 


2006.05.10 


16.9 ks 


501031010 (P 14j 


313.245, 30.859 


59.1 


zUUo.5.12 


lo.Z KS 


^mnQomn /pi ^ 

OUlUo/UlU l^-rlD) 


312.993, 30.653 


59.3 


/UUO.UO.ZO 


17/1 Vo 
i ( .4 KS 


501033010 (P16) 


312.745, 30.450 


59.5 


2006.05.22 


20.0 ks 


501034010 (P17) 


312.207, 30.005 


60.0 


2006.05.22 


13.9 ks 




Lockman Hole 






101002010 


162.938, 57.267 


32.° 7 


2006.05.17 


69.0 ks 



2. Observations and Data Screening 

The observations comprised seven pointings: P8, P12, P13, P14, P15, P16, and P17. 
The FOV of our Suzaku observations are shown in figure 1. We employed revision 1.2 of the 
cleaned event data. We excluded data taken in the low cut-off rigidity < 6 GV. As a background, 
we used a spectrum obtained from the Lockman Hole because its observation date, 2006 May 
17 was close to those of the Cygnus Loop. Obs. IDs, the nominal point, roll angle, observation 
date, and the effective exposure times after the screening are summarized in table 1. Figure 2 
shows a merged XIS1 (back-illuminated CCD; BI CCD) three-color image. Red, green, and 
blue colors correspond to narrow energy bands of 0.54-0.59 keV (O VII Ka), 0.88-0.94 keV (Ne 
IX Ka), and 0.69-0.85 keV (Fe L), respectively We can see strong blue color around the center 
portion (i.e., P15 and P16), whereas red and green colors are enhanced in the NE regions (i.e., 
P8 and P12). For spectrum fitting, we used photons in the energy ranges of 0.2-3.0 keV and 
0.4-3.0 keV for XIS1 and XIS0, 2, 3 (front-illuminated CCD; FI CCD), respectively. 

3. Spatially Resolved Spectral Analysis 

We divided the entire FOV into 119 rectangular small cells, indicated in figure 2, such 
that each cell contained 2500-5000 photons for XIS0 to equalize the statistics. The sizes of small 
cells ranged from 1/32 FOV to 1/8 FOV. We then extracted spectra from them. Generally, 
count rates above 2keV were so low that the statistics were too poor in each cell. It is quite 
difficult to constrain the abundance of S. Therefore, by dividing each FOV into the NE part 
and the SW part, we accumulated photons in the energy range above 2keV from the NE or SW 
half, where each small cell was included, whereas we extracted spectra below 2 keV from each 
small cell. Example regions in each FOV, a black cell for <2keV and a green cell for >2keV, 
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Fig. 2. Three-color image merged by seven XIS FOV (Red: O VII Ka, Green: Ne IX Ka, Blue: Fc L ). 
The data were binned by 8 pixels and smoothed by a Gaussian kernel of a = 25". The effects of exposure, 
vignetting, and contamination are corrected. Small rectangles were the cells where we extracted spectra. 
We show example spectra from black cells for <2keV and green cells for >2keV in figure 3. 
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are shown in figure 2. In this way, we obtained much better constraints on the S abundance for 
each cell than those determined by fitting entire-energy-band spectra extracted from each small 
cell. The S abundances were different for each small cell due to the spectral differences below 
2keV. To generate the response matrix file (RMF) and the ancillary response file (ARF), we 
employed xisrmfgen and xissimarf gen (Ishisaki et al. 2007) (version 2006-10-26), respectively. 
The low-energy efficiency of the XIS shows degradation caused by contaminants accumulated 
on the optical blocking filter (Koyama et al. 2007), which was taken into account generating of 
the ARF file. 

Firstly, all spectra were fitted by an absorbed non-equilibrium ionization (NEI) model 
with a single component [the wabs; Morrison & McCammon 1983 and the VNEI model (NEI 
version 2.0); e.g., Borkowski et al. 2001 in XSPEC v 11.3.1]. Free parameters are the hydrogen 
column density in cm -2 , iVn; electron temperature in keV, kT e ; the ionization time, r, where 
r is the electron density in cm -3 times the elapsed time in sec after the shock heating; the 
emission measure, EM (EM= / n e n-ndl, where n e and nn are the number densities of electrons 
and hydrogens in cm -3 , respectively and dl is the plasma depth in cm); the abundances of C, N, 
O, Ne, Mg, Si, S, Fe, and Ni. We set abundance of Ni to be equal to that of Fe. The abundances 
of other elements included in the VNEI model (i.e., Ar and Ca) were fixed to the solar values 
(Anders & Grevesse 1989). An example spectrum from the black cell in P16, shown in Figure 2, 
is presented in figure 3, indicated as P16 (left) together with the best-fit model. This model 
gave us quite good fits, except for the energy bands around Si- and S-K lines where we can see 
a large discrepancy between our data and the model. Since the emission from both the ejecta 
and the swept-up matter were detected as a projection effect around the center portion of the 
Cygnus Loop (Miyata et al. 1998; Tsunemi et al. 2007; Katsuda & Tsunemi 2008), it is natural 
to consider that we need at least two (i.e., the swept-up matter and the ejecta) components to 
reproduce our data. 

We then applied a two-component NEI model for all spectra. In this model, kT e , r, and 
EM are free parameters for both components. Nr is also left as a free parameter, but tied in the 
two components. Assuming that the swept-up matter surrounds the ejecta, we fixed the metal 
abundances for the swept-up matter component to those of the NE rim regions of the Cygnus 
Loop, where we expect no ejecta component (Uchida et al. 2006). The metal abundances are 
described in the footnote of table 2. Then, we left the metal abundances of 0(=C=N), Ne, Mg, 
Si, S, and Fe(=Ni) for the ejecta component as free parameters. We confined the values of Nn in 
the range from 0.01 to 0.06 xl0 22 cm~ 2 (Inoue et al. 1980; Miyata et al. 2007). Figure 3 shows 
example spectra from seven black cells shown in figure 2. We summarize the best-fit parameters 
for the example spectra in table 2. Based on table 2, we found that the two components clearly 
showed different temperature with each other. The low-/cT e component was responsible for the 
swept-up matter while the high-A;T e component was responsible for another component. We 
found that the metal abundances for at least one element of the high-kT c component were about 
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an order-of- magnitude higher than those for the low-£;T e component, which clearly showed that 
the high-/cT c component represented the ejecta component. In this way, we performed spectral 
fittings for all spectra from 119 small cells. Applying the F-test with a significance level of 
99%, we found that this model gave us better fits than those for one-component NEI models 
for almost all of the spectra [e.g., from figure 3 P16 (left) to P16 (right)]. We obtained fairly 
good fits for all spectra by the two-component NEI model (reduced % 2 < 1.7). If we consider 
the calibration uncertainty of the energy scale (±5eV; Koyama et al. 2007), the values of 
reduced \ 2 become lower than 1.4. Nonetheless, the fits are not acceptable for many spectra 
from a statistical point of view, which suggests that our model is too simple. However, we 
believe that the two-component model is a good approximation to represent our data. Taking 
into consideration that the spectra from almost all cells require the two-component model, the 
ejecta turned out to be distributed inside a large fraction (at least our FOV) of the Cygnus 
Loop. Since the outer edge of our FOV is located around 0.85 R s , our result matches the view 
from XMM-Newton observations that the ejecta are distributed inside ~ 0.85 R s of the Cygnus 
Loop. 

Figure 4 shows maps of the best-fit parameters obtained from the two-component NEI 
model. The subscript, H, denotes the high-/cT e , i.e., the ejecta component while L denotes the 
low-fcr e , i.e., the swept-up matter component. We showed EMs of O, Ne, Mg, Si, S, and Fe for 
the ejecta component. We marked the geometric center of the Cygnus Loop (RA=20 h 51 m 19 s , 
DEC=31°02'48" [J2000]: Ku et al. 1984) as a white dot in each figure. There is a point source, 
1RXH J205036.4 + 302448, within P16 (see, figure 1). We excluded the responsible cell from 
our analysis, which can be seen as a black box in the map of log(rn) shown in figure 4. A^s 
tend to show relatively large values around the center in each FOV to those in the edge regions. 
Since such A^h variations must be due to the contamination problem of the XIS (Koyama et 
al. 2007), the current model of a spatial variation of the contaminants accumulated on the 
optical blocking filter should be modified properly. We found that the value of kT e for the 
swept-up matter component is ~0.2keV in our entire FOV other than P15 and P16, where 
we can see slightly higher values (~0.3keV). The values of kT e for the ejecta component show 
a significant variation from the SW (~0.35keV) to the NE (~0.7keV) portion in our FOV 
as Miyata et al. (2000) already noted the hard X-ray-emitting region in the NE portion by 
ASCA GIS observations and suggested that the hard X-ray might come from the ejecta. The 
collisional ionization equilibrium (CIE) has already been established for both components in 
most cells in PI 7, which seems consistent with the results from Chandra observations of the 
SW rim of the Cygnus Loop (Leahy 2004). 

We noticed a clear unti-correlation between /cT eH and th, and a correlation between th 
and EMs of Si, S, and Fe in the ejecta component. We examined the systematic effect between 
the ionization age (and/or temperature) and the distributions of the Si-, S-, and Fe-ejecta. We 
re-fitted the example spectrum in P15, employing the same two-component model used in this 
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Fig. 3. Top-left panel: X-ray spectra extracted fngn the black cell in P16 indicated in figure 2. The 
best-fit curves for one-component NEI model are shown with solid lines and the lower panels show the 
residuals. The small gap around 2keV is due to the fact that we extracted spectra below 2keV from 
each small cell, while above 2keV from half of each FOV. Top-right panel: Same as top left but with 
two-component NEI best-fit models. The contribution of each component is shown by the dotted lines 



Table 2. Spectral-fit parameters for the example spectra. 



Parameter 


P8 


P12 


P13 


P14 


P15 


P16 


P17 


/VTJx10 22 rm~ 2 l 

iV}] AlU LIll 


n 044+ - 005 

u.U44_ 003 


040+ 005 
u.u^±a_ 003 


0.060±0.004 


0.060±0.004 


49+ - 006 
U - U4y -0.003 


n oro4-o oo** 


n n^i 4-0 nn9 






High temperature component 








/CJ c l[KCVJ 


u.ooitu.ui 


~ rro + 03 


0.60±0.04 


0.60±0.04 


n q o_i_n ni 


n ^q-i-o m 


n Qr;_|_n 01 


U^=L^=1N ) 


u.y±u.i 


U.0±U.l 


47+0-09 


1.2±0.4 


U.oUztU.Uo 


n rc_i_n 1 
U.O±U.l 


U.zU±U.Uo 


Np 


84-0 1 

\J . O— \—\J . -L 


0.7±0.1 


0.8±0.2 


0.7±0.2 


9Q+0 Ofi 


46+ 07 

u - w -0.08 

n 1 4-i-n no 


1 c>+0 05 


A/To- 

lvi g 


U.OZCU. 1 


0.5±0.1 


0.8±0.2 


0.6±0.2 


n i R-r-o no 


n 1 44-0 OQ 


Si 


0.6±0.2 


1.0±0.2 


5.1±0.6 


6.1±0.7 


3.6±0.4 


4.1±0.4 


9 7+0.5 
z - ' -0.6 

44-"? 


G 

o 




2.8±0.5 


10±1 


11±1 


7 1 J_f| Q 




Fpf— Ni 


n fi n+0.06 
U - D/ -0.08 


l.liO.l 


3.3±0.1 


3.9±0.2 




9 Qn-i-n 

Z.OUinu.UtJ 


n 7^,4-0 0*} 


l°g( r l/ cm 3s ) • • ■ 


i i 9 +0.04 
ii ' z -0.09 


1U - ' d -0.04 


11 n +°- 07 
1J - 1J --o.ii 


11.5±0.1 


11 7R+0- 10 


11 7fi+0- 08 


19 0+ 1 - 7 


FlM'f fy 1 D 19 rtn" 5 l 


ri 99-1-0 01 


72+ 006 

u - u,z -0.005 


0.043±0.001 


0.061±0.002 


1 RR4-0 00^! 


n i qi _i_n 009 


o i 094-0 00"} 








Low temperature component 








fcT o2 [kcV] 


0.24±0.01 


0.23±0.01 


0.23±0.01 


0.21±0.01 


0.18±0.01 


0.15±0.01 


0.22±0.01 


Abundances 




(fixed to those determined for the NE 


rim of the Cy§ 


mus Loop)-'- 




log(r 2 /cm~ 3 s) . . . 


11.33±0.04 


11.3±0.1 


11.21±0.05 


11 67+ 09 

-0.08 


13.1±S;? 


12 +1 - 7 


11.58±0.04 


EM 2 [xl0 19 cm~ 5 ] 


4 4+0.2 


9 S+ - 1 


1.91±0.06 


2 16+ 010 

z - iu -0.09 


0.67±0.05 


1 5 +0 - 2 


0.63±0.01 


X 2 /d-oi 


772/643 


749/617 


524/492 


653/576 


654/459 


751/525 


769/523 





* Other elements are fixed to those of solar values. 
The values of abundances are multiples of solar value. 

The errors are in the range A\ 2 < 2.7 on one parameter. 

t r 
EM denotes the emission measure, / n c n^dl. 

*C=0.27, N=0.10, 0=0.11, Ne=0.21, Mg=0.17, Si=0.34, S=0.17, and Fe(=Ni)=0.20 (Uchida et al. 2006). 



paper, but with a fixed ionization age of 1 x 10 11 cm _3 s for the ejecta component (which is the 
typical value in P8, 12, 13). We noticed a slight increase of kT e n that is still significantly lower 
than that in the northern part of our FOV. We found that the emission measures of Si, S, and 
Fe were, respectively, reduced by 30%, 30%, and 15% from those obtained in the case that we 
left the ionization age as a free parameter. These values are still higher than those obtained 
in P8, 12, and 13 by ~50%. Therefore, we can safely conclude that the EMs of Si-, S-, and 
Fe-ejecta in the southern part of our FOV are higher than those in the northern part of our 
FOV. 

The ionization states for the ejecta component in P15 and P16 also reached the CIE 
condition. Apart from P15, P16, and P17, the ionization states for both components are in 
NEI condition in almost all of the cells. We found that all of the elements were distributed 
asymmetric to the geometric center. The EMs for O, Ne, and Mg are enhanced in P8 while no 
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Table 3. Summary of EI (= J n e n,idV for the ejecta component in units of 10 52 cm 3 ) in each FOV.* 



Elements P8 P12 P13 P14 P15 P16 P17 






c Q+0.8 
°-°-0.7 


2.4±0.4 


2.6±0.4 


9 n+ 2 - 3 

z - u -1.0 


2.9±0.6 


3.0±0.5 


1.8±0.2 


Ne 


1.2±0.1 


u.ou_ 08 


0.26±0.08 


0.26±0.12 


0.44±0.09 


0.51±0.08 


0.30±0.05 


Mg 


0.26±0.05 


0.12±0.02 


0.10±0.03 


06+ 05 

u - UD -0.03 


036+ - 038 

U.UOD_ 


0.03±0.02 


0.06±0.03 


Si 


0.25±0.10 


0.24±0.03 


0.49±0.06 


0.7±0.1 


1.3±0.2 


1.2±0.1 


0.5±0.2 


S 


0.18±0.07 


0.23±0.03 


0.40±0.04 


0.63±0.06 


0.9±0.1 


0.9±0.1 


n 1Q+0.10 


Fe 


0.24±0.04 


0.29±0.02 


0.44±0.01 


0.57±0.02 


0.89±0.02 


0.82±0.02 


0.30±0.01 



* Errors quoted are the mean 90% confidence level for each cell. 

enhancements are seen for the other elements there. Also, the EMs for all the elements other 
than Mg show enhancements in P15 and P16. 

We obtained EMs of O, Ne, Mg, Si, S, and Fe for the ejecta component in 119 cells. 
Multiplying the EM by the area of each cell, we obtained the emission integral (hereafter EI, 
EI= / n c ni{dV, dV is the X-ray-emitting volume) of all elements for each cell, and summed up 
the Els for all cells within each FOV. The summed- up Els are summarized in table 3. 

4. Discussion and Conclusions 

We observed the Cygnus Loop from NE to SW with Suzaku in seven pointings. Dividing 
the entire FOV into 119 cells, we extracted spectra from all of the cells, and performed spectral 
analyses for them. Almost all of the spectra were significantly better fitted by a two-component 
NEI model, rather than a one-component NEI model. Judging from the abundances, the high- 
kT e component must be ejecta, while the \ow-kT e component comes from the swept-up matter. 

4-1- Swept-up- Matter Distribution 

The temperature for the swept-up matter component is significantly lower than that for 
the ejecta component. On the other hand, it is similar to that obtained for the rim of the 
Loop, where we expect no contamination of the ejecta (e.g., Miyata et al. 2007). Therefore, we 
believe that we surely separated the X-ray emission of the ejecta inside the Loop from that of 
the surrounding matter. The EM distribution of the swept-up matter (EMl) is inhomogeneous 
in our FOV, as shown in figure 4. The shell of the swept-up matter seems to be thin in the 
SW part (i.e., P15 and P16) of the Loop relative to that in the NE part. Such a trend is also 
reported by XMM-Newton observations covering just north of our Suzaku path (Tsunemi et 
al. 2007). As Tsunemi et al. (2007) mentioned, there might be a blowout in the direction of 
our line of sight around P15 and P16, such as a south blowout of the Loop (see, figure 1). 
The relatively high temperature in P15 and P16 suggests that the velocity of the blast wave is 
higher than that in the other regions. This fact indicates that the density of ambient matter 
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Fig. 4. Maps of the best-fit parameters. The values of TVh, kT e , EMl arc in units of 10 22 cm~ 2 , keV, and 
10 19 cm -5 , respectively. Lower six maps show EMs of O, Ne, Mg, Si, S, and Fe for the high-fcT c component 
in units of 10 14 cm -5 . Wc adjusted the color code, such that we can see the differences in each figure. 
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in P15 and P16 is lower than that in the other regions, supporting the thin swept-up-matter 
shell. 

4-2. Metal Distribution in Ejecta Component 

We divided our FOV into two parts: the NE part (P8, P12, P13, and P14) and the SW 
part (P15, P16, and P17). There is a gap between P16 and P17, which we did not observe. 
Assuming Els of the ejecta in the gap to be averages of those in P16 and P17, we calculated 
the ratios between Els in the NE part and those in the SW part. They are 0~1.2, Ne~1.3, 
Mg~3.2, Si~0.44, S~0.57, and Fe~0.60. For simplicity, we assumed that the mass ratios of 
those elements are equal to the EI ratios. Regarding the O-Ne group, similar amounts of those 
elements exist in the NE and SW parts. In contrast to the O-Ne group, the other elements 
show non- uniform distributions: Mg is distributed more in the NE part by a factor of ~3, while 
Si, S, and Fe are distributed more in the SW part by a factor of ~2. (Note that the mass ratio 
being proportional to a product of the density times the emitting volume will be lower than 
that of the EI ratio.) 

A natural explanation for the asymmetries is an asymmetry at the time of the SN explo- 
sion of the Cygnus Loop. The degree of north-south (1:2) asymmetry of the innermost ejecta, 
such as Si-, S-, and Fe-ejecta, is consistent with that expected from asymmetric explosions 
resulting from hydrodynamic instabilities, which are described in recent theoretical models of 
SN explosions (e.g., Burrows et al. 2007). However, we should keep in mind that there are some 
other possibilities that can produce the asymmetries of the ejecta. The X-ray emission is sensi- 
tive to both temperature and density, i.e., dense ejecta emit more strongly than do thin ejecta 
and cooled ejecta will not continue to emit X-rays. Furthermore, an asymmetric environment 
that might be created by stellar winds (e.g., Blondin et al. 1996; Dwarkadas 2007) as well as 
hydrodynamic instabilities during the remnant evolution ( Jun & Norman 1995) can make the 
ejecta distribution distorted. 

If the mass ejection were to be anisotropic, it would lead a recoil of a stellar remnant 
due to momentum conservation. Therefore, we expect to see the velocity of the stellar remnant 
directed opposite to the momentum of the gaseous SN ejecta (e.g., Scheck et al. 2006). In fact, 
recent observations of a stellar remnant in Puppis A SNR revealed that its velocity was very 
high (>1000kmsec _1 ) and directed opposite to the momentum of fast-moving O-rich ejecta 
(Winkler & Petre 2006; Hui & Becker 2006). Since the Cygnus Loop is believed to be a 
remnant of a core-collapse SN explosion, the presence of a stellar remnant is strongly expected. 
So far, we have not found any stellar remnant associated with the SN explosion of the Cygnus 
Loop (e.g., Miyata et al. 2001). Taking into consideration that the momentum of innermost 
ejecta (which is considered to be strongly related to the stellar remnant) seems to be directed 
toward south of the Loop, we might find the stellar remnant associated with the Cygnus Loop 
in the north of the Loop. We observed the Cygnus Loop from the NE rim to the SW rim 
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with Suzaku as well as XMM-Newton. The observation paths covered about one seventh of 
the entire Loop. Further, Suzaku and/or XMM-Newton observations of the rest of the (11 Oft 3X6 
required to obtain the relative abundances for the total ejecta as well as to reveal the ejecta 
structure in the entire Loop. 

This work is partly supported by a Grant-in- Aid for Scientific Research by the Ministry 
of Education, Culture, Sports, Science and Technology (16002004). This study is also carried 
out as part of the 21st Century COE Program, 'Towards a new basic science: depth and 
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